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THE STRESS-STRAIN BEHAVIOR OF POLYMER-NANO TUBE COMPOSITES 
FROM MOLECULAR DYNAMICS SIMULATION 


S. J. V. FRANKLAND 1 ’ 2 , V. M. HARIK 1 ’ 2 ’ 3 , G. M. ODEGARD 1 ’ 2 , D. W. BRENNER 4 , AND T. S. GATES 5 

Abstract. Stress-strain curves of polymer-carbon nanotube composites are derived from 
molecular dynamics simulations of a single-walled carbon nanotube embedded in polyethylene. A 
comparison is made between the response to mechanical loading of a composite with a long, continuous 
nanotube (replicated via periodic boundary conditions) and the response of a composite with a short, 
discontinuous nanotube. Both composites are mechanically loaded in the direction of and transverse to 
the NT axis. The long-nanotube composite shows an increase in the stiffness relative to the polymer and 
behaves anisotropically under the different loading conditions. The short-nanotube composite shows no 
enhancement relative to the polymer, most probably because of its low aspect ratio. The stress-strain 
curves are compared with rule-of-mixtures predictions. 
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7 1NT - Elastic- modulus of the nanotube 

Y 2 - Transverse elastic modulus of composite 

7 2 nt - Effective transverse elastic modulus of nanotube 

Y m - Elastic modulus of the polymer matrix 

8 y - Strain tensor 

Sjj - Potential well depth 

(j>(d) - Torsional potential for angle 6 

<1 ) a (r) - Potential energy at the atom a for separation distance r 

v a - Magnitude of the velocity of atom a 

v“ - /-component of the velocity of atom a 

6 - Torsional angle around the CH 2 -CH 2 bond 

a y - Stress tensor 

cf u - Van der Waals separation distance 


1. Introduction. In the last few years, single-walled carbon nanotube-polymer composites 
have generated considerable interest in the materials research community because of their potential for 
large increases in strength and stiffness when compared to typical carbon-fiber-reinforced polymer 
composites. Even though some nanotube composite materials have been characterized experimentally [1- 
12], the development of these materials can be greatly facilitated by using computational methods that 
allow for parametric studies of the influence of material and geometry. Molecular dynamics (MD) 
simulations can predict the effect of mechanical loadmg on specific regions of polymer-nanotube (NT) 
composites. In this work, stress-strain curves are generated from molecular dynamics simulations of 
r\r\l ymer-NT composites to explore the nano-structural effects on the overall mechanical properties of 
polymer-NT composites. 

Various aspects of the mechanical reinforcement of polymers by carbon nanotubes (NT) have 
been addressed computationally. One example of recent work is an equivalent-continuum modeling 
method that was used to predict the elastic properties of two single-walled polyimide-NT composites for 
various NT volume fractions, lengths, and orientations [13]. At the continuum level, finite element 
analysis has been used to predict the macroscale effects of ‘wavy’ or not straight NTs [14, 15], Molecular 
statics simulations of the transverse mechanical properties of NT bundles have shown the strong cohesion 
of NTs compared to a graphite-epoxy matrix [16]. Mechanical properties of NT yams have been 
predicted from the elastic constants of NT bundles [17]. Molecular modeling has predicted the adhesion 
between a NT and various polymer matrices [18]. Molecular dynamics simulations have predicted that 
the polymer-NT interface can be reinforced by covalently bonding the NT to the polymer matrix [19]. 

Another feature of the mechanical behavior of polymer-NT composites that can be modeled is the 
stress-strain relation. The stress-strain relationship provides the overall mechanical response of NT- 
polymer composites when subjected to mechanical loading. In the present work, stress-strain curves of 
two polyethyl ene-NT composites are generated with molecular dynamics (MD) simulations [20]. The 
effects of both the anisotropy of the composites when subjected to different loading conditions and of the 
length of the NT are explored. The first composite contains a long (infinitely long) NT, and the second 
composite, a short NT. Details of the composite structure and MD simulations are given m Section 2. 
Stress-strain curves are presented for longitudinal and transverse loading conditions of each composite 
(Section 5). A detailed description of the computation of stress from molecular force fields in Section 3. 
Comparison is made of the simulated stress-strain curves with the rule of mixtures. The relevant 
equations are in Section 4. 
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2. Molecular Dynamics Simulations. The unidirectional NT-polymer composites contain 
symmetrically placed long or short NTs as shown in Figure 1. The dashed boxes in Figure 1 enclose a 
representative volume element that is simulated by MD. In each composite, the NTs are sufficiently 
separated by polymer to prevent direct NT-NT interactions. The NT composite in Figure 1(a) contains a 
periodically replicated (10,10) NT which spans the length of the simulation cell. (The (n,m) notation 
refers to the chiral vector of the NT in terms of the primitive in-plane lattice vectors of a graphene sheet 
[21].) In this composite, the NT is embedded in an amorphous polyethylene (PE) matrix, which is 
represented by beads of united atom -CH 2 - units. Specifically, the PE matrix has eight chains of 1095 - 
CH 2 - units. The short-fiber composite, shown in Figure 1(b), contains a 6-nm capped (10,10) NT which 
is approximately half the length of the simulation cell. The NT caps each consist of one half of a C 240 
molecule. In this composite, the amorphous PE matrix contains eight chains of 1420 -CH 2 - units. The 
overall dimensions for the unit-cell, or minimum image, of each composite in the MD simulation are 
approximately 5x5x10 nm. Periodic boundary conditions were used to replicate the cell in all three 
dimensions. For comparison to the NT/PE composites, an equivalent-sized block of amorphous PE 
without a NT was also simulated. For each structure, a PE density of 0.71 g/cm 3 was used. 

In the MD simulation, the van der Waals interfacial interaction between the polymer and the NT 
was modeled with the Lennard-Jones potential [22], 

frrl 2 * r 6 ^ 

(2.1) U = 4s t 
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where U is the potential energy between a pair of atoms, r is the separation distance between the pair of 
atoms, £u is the potential well depth, and Ou is the van der Waals separation distance. For the interaction 
between the carbon atoms of the NT and the polymer -CH 2 - units, the potential was parametrized with 
^ ,= 0.4492 IJ/mol and 0.3825 nm [22,23]. The PE chains were simulated with a molecular- 
mechanics force field adapted from Reference 23. Specifically, the -CH 2 - units of the PE chains were 
separated by bond lengths of 0.153 nm by using the SHAKE algorithm, a constraint dynamics method 
which constrains the bond length within a user-defined tolerance [22]. Angle bending forces were 
modeled with a harmonic valence-angle potential having an equilibrium angle of 1 12.813° and a barrier of 
520 kJ/mol. A torsional potential <fi(0) of the form 


(2.2) <f>(0) = Ai( 1 + cos Q) + A 2 (1 + cos 20) + (1 4- cos W) 

was used where 6 is the torsional angle around the CH 2 -CH 2 bond, and Aj— 1 8.087, A 2 =-4.88, and 
A 3 =31.8 kJ/mol [23]. The Lennard-Jones potential was also used to describe non-bonding interactions 
between -CH2- units in either the same chain or between different chains. For these interactions, 
8ij=0.4742 kJ/mol and Ou=0.428 nm [23]. The NT was modeled with a many-body bond-order potential 
developed for carbon [24]. This carbon potential is parametrized for C-C bonds of lengths up to 0.17 nm 
which is within the magnitude of the strain applied to the composites in the present work. 

The MD simulations were carried out by using DL-POLY [25], a large-scale MD simulation 
package available from Daresbury Laboratory. This program was adapted to include the carbon potential 
for the NT and to simulate the application of strain to the composite. All simulations were carried out at 
300 K, with a 2 fs timestep. 


3. Stress-strain Curves from Simulation. Stress-strain curves were generated for the long- and 

short-NT composites and for the pure polymer via MD simulation. For both composite configurations, the 

longitudinal (parallel to the NT axis) and transverse responses were simulated. 


3.1. Strain. For each increment of applied deformation, a uniform strain was prescribed on the 

entire MD model. For the longitudinal and transverse deformations, pure states of the strains an ^ s 2 2 , 

respectively, were initially applied (see Figure 2). The application of strain was accomplished by 
uniformly expanding the dimensions of the MD cell in the direction of the deformation and re-scaling the 
new coordinates of the atoms to fit within the new dimensions. After this initial deformation, the MD 
simulation was continued and the atoms were allowed to equilibrate within the new MD cell dimensions. 
This process was carried out for the subsequent increments of deformation. The applied strain increment, 
in either the longitudinal or transverse direction, was 2%, and was applied in two increments of 1% each 1 
ps (500 steps) apart. After each 2% increment of strain, the system was relaxed for 2 ps, then the stress 
on the system was averaged over an interval of 10 ps. For each configuration, six increments of 2 % 
strain were applied up to a total of approximately 12 % over a period of 72 ps. The corresponding strain 
rate was l.OxlO 10 s' 1 . This high strain rate is inherent to MD simulation which includes dynamical 
information usually on ps or ns timescales. 

3.2. Stress. In general, the stress in a solid (or a group of interacting particles in the form of a 
solid) is defined as the change in the internal energy (in the thermodynamic sense) wdth respect to the 
strain per unit volume. For example, at the continuum level, the stress tensor for a linear-elastic material 
is [26] 


(3.1) 
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where V is the volume of the solid, E is the total internal energy, s^- is the strain tensor, and the subscript S 
denotes constant entropy. Assuming that the internal energy is equal to the strain energy of the solid, then 
Hooke’s law may be derived from equation (3.1). Furthermore, if the strain energy is expressed in terms 
of an allied force acting over the surface area of a solid, then a more familiar form of stress as force per 


unit area is derived. 

At the atomic level, the total internal energy given in equation (3.1) can be expressed as the 
summation of the energies of the individual atoms, that compose the solid: 


(3.2) 


E a = T a + U a 


i Af“(v“) 2 + ®“(r) 


where for each atom a, T* is the kinetic energy, If 1 is the potential energy, Af" is the mass, v a is the 
magnitude of its velocity, and O a (r) is the potential energy at the atom location r. Using a Hamiltonian 
based on these individual energy contributions, E a , it has been shown that the stress contribution for a 
given atom is 


(3.3) 


<*=-■ 


i 

U 


( ' 

1 + Y J F, a *rf 

S. P 

where V 1 is the atomic volume of atom a, v“ is the /-component of the velocity of atom a, y a , is the j- 
component of the velocity of atom a, ir ap is the /-component of the force between atoms a and P 
obtainable from the derivative of the potential <X>(r), and rj® is the /-component of the separation of 
atoms a and p (Figure 3) [27,28]. 

The stresses that were used to generate the stress-strain curves for the NT composite were 
average atomic stresses for the volume of the model. Therefore, the stress components of each model 
were calculated for each strain increment by using 


(3.4) 
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where V is the volume of the MD model and V = i V a . The stress calculated with equation (3.4) was 
then averaged over time via the MD simulation. ° 

4. Rule-of-Mixture Analysis. For a nanocomposite under uniaxial loading, the dependence of the 
elastic modulus on the NT volume fraction can be estimated by the rule of mixtures [29]. The 
longitudinal elastic modulus, Y\, of the composite cell with long NTs under isostrain conditions is 


(4.1) 


r, = W«- 


+ YJ m 


where Tint and T m are the effective longitudinal elastic modulus of the NT and the polymer matrix, 
respectively, and fm and f m are the volume fractions occupied by the NT and the polymer matrix, 
respectively. The transverse modulus of the nanocomposite, T 2 , is 


(4.2) 
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where F 2 nt is the effective elastic modulus of a NT in the transverse direction. The volume fractions 
satisfy 


(4.3) 


fm + fm - 1 - 


The (10,10) NT used in the present work has a radius small enough to be thought of as a solid beam 
[30]. Therefore, its volume fraction, includes the entire NT cross-section and is given by 
7T a. „ Yx/Vi^rp' Vi .... io tUp pnmliKrinm tron rlpr \X7Valc ept^aT*ati£vn dlStanCC betWCCU the NIT and 
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the matrix, and A ce n is the cross-sectional area of the unit cell. The van der Waals separation distance 
depends on the nature of the polymer/NT interfacial interactions. 


5. Results and Discussion. The stress-strain curves of the long- and short-NT composites 
generated from MD simulation are presented in Figures 4, 5, and 6. Before describing the individual 
stress-strain curves in more detail, some general comments are in order. First, the endpoints of the curves 
are indicative of how far the simulation was run, and hence, should not be taken as yield points. No 
attempt is made in this work to simulate the yielding of carbon nanotubes. Also, for both composite 
configurations, the maximum strains are on the order of 10 % after the unit cell had been subjected to 
stepwise increases of applied strain. The strain leads to an overall displacement of 0.5 nm for each unit- 
cell edge in the longitudinal direction and 0.25 nm in the transverse direction. 

In Figure 4, the stress-strain curve (solid triangles on solid line) of the long-NT composite under 
longitudinal loading conditions is compared with the rule of mixtures (solid line without symbols) from 
equation (4.1). Of the composites simulated here, this one has the most significant enhancement in its 
stress-strain curve relative to the polymer (to be shown in Figure 5). For this sample, both the NT and 
the polymer were subjected to constant strain conditions. The stress-strain curve is comparable at low 
strain to the stress-strain behavior predicted by the rule of mixtures in equation (4.1) using a volume 
fraction of 8%, which is the NT volume fraction in the simulated composite, and an effective modulus of 
600 GPa for the (10,10) NT [31]. 

The stress-strain curve of the short-NT (solid circles on solid line) subjected to longitudinal 
loading conditions is plotted with the simulated polymer stress-strain curve (solid squares on solid line) in 
Figure 5. Almost no enhancement relative to the polymer is observed in the stress-strain curve of the 


composite with the short 6-nm NT. The bulk polymer stiffness is improved by less than a factor of two 
when the NTs are added into the polymer. This lack of effect is consistent with the low aspect ratio of the 
nanotube, 1:4, which is considerably lower than the more than 1:1000 aspect ratio expected for a typical 
single-walled NT. Length scales of this magnitude are not accessible with the current MD simulation. 
Furthermore, the simulated interface is established by non-bonding van der Waals forces between the 
polymer and the carbon NT [32], and there are no chemical bonds [19] or other strong interactions [13] to 
strengthen the mterfacial adhesion. 

The stress-strain curves of both composites subjected to transverse loading conditions are 
compared with each other and the polymer (solid circles on solid line) in Figure 6. All three MD 
generated curves are similar. For the composite with the long NT (solid squares on solid line), the stress 
after loading m the transverse direction is approximately 30 times lower than the stress levels after 
loading to the same strain level in the longitudinal direction (Figure 4). This difference in longitudinal 
and transverse behavior emphasizes the anisotropy of the composite. Both the long- and short-NT 
composites have more similar transverse behavior than the polymer which may be because the NT fills 
the same relative cross-sectional area in both composites. The solid line in Figure 6 is from the transverse 
rule of mixtures given by equation (4.2) for transverse loading conditions of the long-NT composite. For 
this comparison, an effective NT transverse modulus of 10 GPa [33] is used, and a polymer modulus of 
2.7 GPa, which is derived from the initial points of the simulated polymer stress-strain curve. For 
somewhat higher estimates of the NT transverse modulus, the predictions will not change because the 
contribution from the matrix modulus is dominant in equation (4.2). 


6. Concluding Remarks. Molecular dynamics (MD) simulations were used to establish the 
stress-strain behavior of polymer-carbon nanotube (NT) composites mechanically loaded in both the 
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discontinuous fibers were investigated. Stress at the atomic level was defined by using energy 
considerations, averaged over the volume of the MD model. The simulated stress-strain curves were 
compared with the rule of mixtures. 

The stiffest behavior, and the one closest to approaching the rule of mixtures, was observed for 
the longitudinally loaded, long continuous NT composite. This composite was very anisotropic shoving 
stresses that, when subjected to transverse loading conditions, were more than one order of magnitude 
lower than when subjected to longitudinal loading conditions. By contrast, the short, discontinuous fiber 
composite, upon loading in either direction, showed no appreciable load transfer from the polymer to the 
NT. It is anticipated that increasing the aspect ratio of NT in this composite would lead to further 
enhancement of the stress-strain curve and greater differences between the longitudinal and transverse 
stress-strain behavior [19,34]. 
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FIG. 1 . Schematic of polymer nanocomposites filled with long and short carbon nanotubes 







FIG. 4. Longitudinal stress-strain relation oflong-NT composite and comparison to rule of mixtures 
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FIG. 5. Longitudinal stress-strain relation for short NT composite compared with the stress-strain behavior of the polymer 
matrix and the rule of mixtures for transverse modulus 
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